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Noise in optical Telecom fibers is an important limitation on optical quantum data transmission. 
Unfortunately, the classically successful amplifiers (such as EDFA) cannot be used in quantum 
communication because of the no-cloning theorem. We propose a simple method to reduce quantum 
noise: the insertion of phase-shifters and/or beam-splitters at regular distance intervals into a fiber. 
We analyze in detail the case of qubits encoded into polarization states of low-intensity light, which 
is of central importance to various quantum information tasks, such as quantum cryptography 
and communication. We discuss the experimental feasibility of our scheme and propose a simple 
, experiment to test our method. 
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I. INTRODUCTION 



Quantum communication (QC) has recently emerged as a subject of much interest, due to its applications in dis- 
tributed quantum computation and quantum cryptography l]. In QC, non-orthogonal quantum signals are typically 
I i transmitted through Telecom fibers. Reducing noise in Telecom fibers is crucial for QC applications, because the 
I '■ very weak signals carried by polarization states are usually employed. Ideally, a single photon in a superposition of 
two pulses separated in time with a controlled phase difference (i.e., |pulse A) -\- e*^|pulse B)) may be used. How- 
PsJ ' ever, on-demand single-photon sources remain an important technological challenge. Currently, weak coherent states 
K*" ', are often employed as approximate single photon signals. It is well known in quantum cryptography 0i 0) ~ ^ 
OO branch of QC - that weak coherent states may open up loopholes in security because of the probabilistic existence 
' of multi-photon signals Indeed, a weak coherent state, when phase randomized via decoherence, gives a Poisson 
I distribution in photon numbers. An eavesdropper. Eve, may, in principle, measure the photon number in such a signal 
' by a non-demolition measurement. Afterward, she can stop all single photon signals from reaching the receiver. Bob. 
For each multi-photon signal, she can steal one photon and keep it in her quantum memory and send the rest of the 
signals to Bob by using, for example, an ideal channel instead. Since Eve now has an exact copy of the quantum state 
transmitted to Bob, this creates a significant challenge in the security of quantum cryptography 0. Thus, attenuation 
^ ^ losses and decoherence in QC are a major issue and methods for reducing such quantum noise are therefore important. 
I \ Unfortunately, the classically successful amplifiers (such as EDFA 0) cannot be used in QC because of the no-cloning 
^ . theorem @, and new methods must be explored. 

d ' Here, we propose a method to reduce noise in the transmission of quantum optical signals in a Telecom fiber. Our 
^ method is inspired by the theory of quantum dynamical "bang-bang" (BB) decoupling 7J. However, a key novelty of 
our work is the following: we propose to implement BB control in space, rather than time, through the insertion at 
^ , regular intervals of a sequence of simple linear optical elements (phase-shifters and/or beam-splitters) in sections of 
a Telecom fiber. We also discuss the experimental feasibility of our scheme, and propose a few simple experimental 
. ^ tests. We do not expect our method to improve the fidelity of classical light transmission compared to, e.g., EDFA 
' amplifiers, since our method turns out to be quite sensitive to reflection from optical elements and deviations from 
" " ' average fiber homogeneity, which is not the case for classical amplifiers. 



II. QUANTUM NOISE IN OPTICAL FIBERS 



An optical fiber provides boundary conditions that guide light along a (locally) straight trajectory. An ideal 
fiber allows modes of traveling photons to propagate through unchanged. A real fiber induces noise (dispersion, 
loss, decoherence) compared to the ideal case. The method we introduce in this paper is designed to cancel (in 
principle) all quantum noise. The dominant classical-light loss mechanisms in an optical fiber are UV absorption, 
Rayleigh backscattering, OH absorption, and Raman scattering, and lead to typical attenuation rates, for state-of- 
the-art commercial silica Telecom fiber, of about 0.25db/km |5(. These mechanisms are active also in the quantum 
regime 0, 0| . All noise processes affecting quantum light in optical fibers are derivable from a microscopic Hamiltonian 
describing (i) the direct interaction between photons and the optical (dielectric) material of a fiber, and (ii) the indirect 
interaction between photons and quasi-particle excitations of the fiber material, such as polaritons and photon-phonon 
coupling. These indirect interactions are, of course, in turn derivable from a microscopic Hamiltonian that takes into 
account matter-matter interactions in the fiber, and couples them to photons. The derivation of the resulting effective 
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interactions (e.g ., a non-linear Schrodinger equation) from such microscopic Hamiltonians has been covered in detail, 

e.g., in liiEa- , 

The starting point of our analysis is the observation that all interactions involving photons can be written in terms 
of polynomials in the bosonic raising and lowering operators bj , bj (where j is the mode of the traveling photons 
0, 0]). A polynomial of order N describes an iV-photon process, and typically the cross-section of interactions 
decreases with increasing N. In the case of the non-relativistic quantum electrodynamics of charged particles one can 
decompose the photon-matter interaction Hamiltonian into linear and quadratic terms with respect to the photon 
field, Hi ~ H\ + Hf, where the linear part is 

H\ = Y,{b,B]+b]B,), (1) 

j 

where the "bath" operators Bj depend only on the variables of charged particles and/or quasi-particles, and the 
quadratic part is a function of the bilinear operators, b\bj, and bibj. Higher-order interactions may originate 
from relativistic effects. In general iJj, which makes no contribution to one-photon processes, is much smaller than 
H\ Therefore, the quadratic term can usually be neglected. 

Let us substantiate these arguments by briefly reviewing the corresponding non-relativistic electrodynamics. Con- 
sider particles a with charge qa and mass ttIq, which constitute the optical material of a fiber. Let Tq, and Pa be 
the position and momentum of particle a, and A(r) be the vector potential of the photon field. The system-bath 
Hamiltonian that describes the dynamics can be written, in the Coulomb gauge, as 

H = Hf) + Hi; 

Hf) = Hm + Hp (2) 

Here Hm depends only on the variables of the charged particles. Hp ~ + 1/2) is the free photon 

Hamiltonian, where bj, 5^ are the photon annihilation and creation operators in the normal vibrational mode j of the 
field identified by the wave vector kj , the polarization Sj and the frequency ujj = ckj , where c is the speed of light in 
vacuum. Then the linear part with respect to the photon field Hrj is 



H\ = P" • A(ra) + f^S„ . B(r„)) (3) 

nia Inia 



]{B]b,+B,b]) 

where for a cubic box with dimension L the operator Bj can be expressed as 



j 



which only depends on the variables of charged particles. Here ga is the g factor, Eq is the permittivity of free space, 
and Sq is the spin of particle a. Note that the interaction is linear in the operators bj and 6j. 
The quadratic part of the interaction Hamiltonian is found to be 

a 

and is a function of the bilinear operators, blbj, b^b'j and bibj . 

Under the long-wavelength approximation, where the spatial variations of the electromagnetic field over the size 
of the particles is negligible, first order perturbation theory of ijj results in the widely applied dipole interaction 
(e.g., jlOi | and references therein). Some effective interactions, such as atom- mediated photon-photon interactions 
and nonlinear photon- photon interactions (Kerr effect) have been derived without consideration of iJ| 0, 0,^3- 
provide more details on these effective interaction in Section [VIII For simplicity of presentation we will first design an 
"anti-Zmear-decoherence fiber" by considering ffj only. Later on we show how to treat higher order interaction terms. 
It is important to stress that in essence our method hardly depends on the details of the interaction, but depends on 
the statistics of photons as bosons. For this reason our method is very general and is in principle applicable to the 
entire phenomenology of quantum noise processes affecting photons in fibers, though its practical applicability is a 
matter of being able to satisfy certain constraints that will be discussed in detail below. 
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III. THE ANTI-LINEAR-DECOHERENCE FIBER 



We first consider quantum data transmission through a Telecom fiber with noise induced by H\. Since H\ describes 
absorption and creation of photons, it generates photon loss, among other processes. To simplify, we suppose that a 
polarization photon is transmitted from end A to end B. One can define a logical qubit supported by |0)^ = b\ |vac) 
and — h\ |vac) where the mode indices refer to the two polarization states. The initial state at end A is 
l^*^) = (a|0)^ + h\l) j^) \M), where \M) is the state of the bath (dielectric material and quasi-particle excitations 
in the fiber). At the time T — X/v (where X is the distance between A and B, and v is the average speed 
of light in the Telecom fiber) the wave function is \'^{T)) = U{T,0) where the evolution operator is (in 

units where h = 1) C/(T,0) « g^^-f^^^'^)^ • • • e-*-^(2A)rg-,//(A)r ^ .^j^^^.^ H{kA) = s /(fc^ijAl^^^;) + Ho{x)]dx is 

the average Hamiltonian over the fcth segment, where Hq is a sum of the matter (and/or excitations) and photon 
self-Hamiltonians, r — A/v, and we have assumed that N = X/A is large in order to expand the normal-ordered 

exact propagator U{T,0) = : exp[— i J^[Hj{x) + HQ{x)]dx] :. I.e., we have neglected deviations from average fiber 
homogeneity, Sk — {{H{kA) — [Hi{kA) + Ho{kA)]y) [t/(T, 0) can easily be expressed including such second and 
higher order moments using a Magnus expansion, and it is known how to generalize BB decoupling to treat such 
higher moments, at the expense of more BB pulses ^2]- The interaction Hj entangles the output wave function 
at end B with the material/excitations in the fiber. By standard arguments it follows that, therefore, the quantum 
information encoded into the photon state will decohere 

In order to solve this problem of decoherence, we draw inspiration from the idea of BB decoupling via time-dependent 
pulses 7] (we note that a method for finding such pulses directly from empirical data was proposed in [T^'l. We first 
recall the action of a phase-shifter. It is simple to show [using the Baker-Campbell-Hausdorff (BCH) formula 0|] for 
a boson that 

gi^n^te-i-^" = e^'^6^ e''^"6e-*'^" = (5) 

where n = b^b is a boson number operator. Physically, the operation e"" is a tt phase-shifter (it puts a phase of tt 
between the number states |0) and |1), not to be confused with our logical qubit states). Defining the 7r-phase-shifter 
operator 

n = = e''^("i+"^\ (6) 

we therefore have 

UHU = Ho-Hl (7) 

because the photons term of Hq is ^ ?iiOj{nj + 1/2), so that [Hq^ui + n2\ = 0. The crucial point is that the sign of 
the linear term of the interaction Hamiltonian has been negated by the action of two phase-shifters, i.e., effectively 
time-reversed. Now, if we install thin phase-shifters inside the fiber at positions x — 0, A,2A, from A to B, the 
evolution will be modified to 

= [iv,...,n,2,n,i,n], 

where in the second line we have introduced a self-explanatory notation that will be used repeatedly below. Note that 
in writing this expression we have neglected the variation of H inside the phase-shifter; this will hold provided that 
the phase-shifter width is much smaller than the distance over which deviations Sk from average fiber homogeneity 
become significant. Further note that we are applying the "parity- kick" version of BB decoupling but are 

implementing it in space, rather than time. Now assume that the average Hamiltonians over two successive segments 
are equal: 

H\{{k + 1)A) = H\[kA) 

Ho{{k + l)A) = HoikA). (8) 

The better this approximation, the better our method will perform; we address deviations in Appendix fXl In this 
case, to first order in r, and using Eq. ((7|), we have exact cancellation of ijj between successive segments: 

^-iH({k+l)A)TYip-tHikA)TYi ^ ^-iH{{k+l)A)T ^-inH(kA)UT ^ ^-2iHo{kA)T j^g-j 

This yields the overall evolution operator 

U'{T,0) = f,-iHo{X)NT ^ ^-iHo(X)T^ 

i.e., the evolution is completely decoherence- free, in analogy to the ideal BB limit of infinitely fast and strong pulses 
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IV. ROUGH ESTIMATE OF REQUIRED INTER-PHASE-SHIFTER DISTANCE 



Because of the in-principle equivalence between the BB method and the quantum Zeno effect [l6j, the proposed 
method can only work if the phase-shifters are inserted at small intervals A over which coherence loss is quadratic 
( "Zeno-like" ) , rather than exponential ( "Markovian" ) . A reliable estimate of A requires a first principles calculation 
which is beyond the scope of the present work; we present a phenomenological model for a detailed estimate of A in 
Appendix ^ Here we give a rough upper bound estimate of this distance. We assume that the linear term of the 
interaction Hamiltonian gives rise to the 0.25dB/km (5 x 10~^) classical loss figure in a Telecom fiber. Our main 
approximation now consists in further assuming that the insertion of phase-shifters into the fiber causes a reduction 
of loss from first to second order, and we use this to estimate the A required in the quantum case. Thus, imagine 
a distributed quantum computing scenario where small-scale quantum computers are connected by optical fibers of 
length about 1km. Our goal is to have reliable quantum computation within the fault-tolerance threshold value of 
10~^ error rate for each elementary quantum logical operation. [We remark that for reliable quantum communication 
of entangled photon pairs, the current error rate of about 5 x 10~^ is already acceptable provided one allows the 
application of entanglement purification p^ : our scheme is significantly simpler.] Therefore, we need to cut down the 
loss figure from 5 x 10~^ to say 10"'*. Suppose we need to insert N phase-shifter within 1km of a Telecom fiber. Denote 
the attenuation between a pair of phase-shifter by I. Then, without the N phase-shifters, we have (1 — l)^ = 0.95. 
For a sufficiently large TV, we can expand the expression binomially and obtain the approximation IN = 0.05. Now, 
with the insertion of phase-shifters, we simply assume that the attenuation between two phase-shifters is due to a 
second order contribution of the form P. We further assume that those contributions sum up in usual addition. 
Therefore, we have I'^N = 10~^. This yields I = 2 x 10~^ and N = 25. Recalling that two phase-shifters are needed 
per cancellation step, we see that about 50 phase-shifters have to be inserted in a distance of 1km which translates 
to one phase-shifter every 20m. This figure is merely a rough upper bound estimate on the distance A between two 
phase-shifters for our scheme to be useful; one can also determine A via the experiment we propose below. Also note 
that we have assumed here that the fiber is straight as is typically done in theoretical models. In order to regain the 
straight fiber approximation, in the case of a curved fiber A is upper-bounded by the local radius of curvature. 

While in spirit our method is similar to BB decoupling 0, a major advantage here is that we do not need to 
apply any time-dependent pulses, which may result in significant uncertainties such as gate errors and off-resonance 
transitions. Instead, the phase-shifters may be incorporated into the fiber directly during the manufacturing process. 
Alternatively, time-independent (say, electronic or pressure) controls may be applied at various points of a Telecom 
fiber to achieve the action of pulse shifters. 



V. THE ANTI-BILINEAR-DECOHERENCE FIBER 



We now consider higher order processes. Although they are generally weak, the bilinear interactions appearing 
in Hj may still cause decoherence. A direct harmful consequence is to change the polarization direction, through a 
term such as b\b2- In the classical case, the fiber structure can be designed so that a known polarization direction 
can be preserved p^ . In the quantum case the polarization direction is not known prior to the transmission and the 
classical method is not applicable. In this case one must in general consider a system-bath Hamiltonian that is a linear 
combination of all 10 possible independent bihnear terms: A = {b\b2, 62^1, {b\)^, (^2)^; (^i)^; (^2)^}, S = {6162, ^i^2}i 
C = {b\bi,blb2} (the grouping will be clarified momentarily). It can be shown that all 10 of these terms can be 
eliminated by installing 18 linear optical devices that include beam-splitters in addition to phase-shifters, i.e., in 16 
elementary steps (we combine beam-splitting and phase-shifting into one step). This result is based on Eq. (jS)) and 
the following identities [that follow directly from Eq. (0] 

e#ft(5t)2g-»0rl ^ e2#(^t)2^ g#n(^)2g-.0« ^ e-^'''' (b)^ . (10) 

The role of the beam-splitter is to eliminate the set of operators C; the beam-splitter is inserted after the first eight 
steps. The 16-step result can be considerably simplified in the realistic situation wherein the two polarizations used 
to represent our qubit are degenerate. In this case C becomes b\bi + 62^2, which generates an overall phase and hence 
will not cause decoherence. In this degenerate case, as we now show, we need only phase-shifters to eliminate all 
contributions to decoherence. Let 

H, ^e*'^"', r = e"^^^-^^'>/^, (11) 
i.e., a pair of phase-shifters. It follows immediately from Eqs. (|3|l. ((Tni) that 
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while 

T^Br = B, nUn = A u^Bn = B (12) 

[where 11 = 111112 was used above]. From these and the results for the "anti-linear-decoherence fiber", the sequence 
ili2 = [2. n, 1, n] does not contain any linear terms, but still contains all bilinear terms. Then, the sequence 

^^l234 = [r!34,rt,r!i2,r] 

= [4,n,3,nr^2,n, i,nr] (13) 

has, in four elementary phase-shifter steps, eliminated Hj as well as A, and in particular the polarization-direction- 
changing terms b\b2 and 62^1: at this point we have a polarization-preserving fiber. Note that the composite terms 
can be combined into a single phase-shifter, i.e., 

Jjpt _ g'i7r(ni+3fi2)/2 
nr — g»'^(3"i+n2)/2^ 

The only remaining bilinear terms at this point are the counter-rotating terms B = {6162, 62}' that are typically 
neglected in the rotating wave approximation To eliminate them nevertheless, we note that 

UiBUi ^ -6. 

Therefore the sequence that eliminates all linear and bilinear terms for degenerate qubit states is 

[1^5678, Hi, r!i234, Hi] = [8, H, 7, nrt, 6, n, 5, nrni, 4, n, 3, nrt, 2, n, 1, nrni], 

which involves 8 elementary phase-shifter steps (note that Xirili = e'^(^"i+"2)/^). At this point we have a fiber that 
is completely free of both linear and bilinear decoherence-causing terms for degenerate polarization qubits. 

We can repeat the mixed-classical-quantum rough distance estimate above, by simply assuming that now contribu- 
tions to decoherence come only due to third order in I: l^N = 10~*. This leads to = 5/\/20 « 1.2, and recalling 
that 8 phase-shifters are needed per cancellation step, we arrive at an upper-bound estimate of about 10 phase-shifters 
per km, or one phase-shifter every 100m. These phase-shifters must be introduced in addition to the ones used above 
for cancellation of first order effects. We have again assumed here that the fiber is straight; local curvature may 
impose a lower upper bound. 

VI. GENERAL DECOHERENCE ELIMINATION 

So far we have considered linear and bilinear photon terms in the interaction Hamiltonian. The most general 
two-mode photon-related term in a Hamiltonian is 6|'^6i62'^&2- Provided r ^ s and k ^ I the identity 

shows that such a term can be eliminated using only phase shifters. For example, when r-|-s-|-fc-f-Hsan odd number, 
our considerations in the linear case show that the term can be eliminated using the phase shifter 11, while b^x^\ can 
be eliminated using e"*^"^. High-order terms with r, s, /c, ^ > 1 arise if one considers the relativistic contribution, and 
they appear also in most of the effective photon scattering theories. It should be clear that if such terms arise, they 
can be reduced using additional phase-shifters, or beam-splitters in the case r = s and/or k = I, which arise due to 
terms containing photon number operators. 

VII. CONNECTION TO KNOWN LEADING LOSS MECHANISMS IN OPTICAL FIBERS 

As mentioned in Section ^] the leading loss mechanisms in optical fibers are well characterized: UV absorption, 
Rayleigh backscattering, OH absorption, and infrared absorption. It is useful to quickly review how these processes 
arise, and then are treated by our method. Consider, for example, the case of Rayleigh backscattering. We base 
our discussion on the standard reference (for a general description of absorption see p. 168; the cross section of 
Raleigh scattering is given on pp. 371-373). The discussion starts 19] [Eq. (4.9.9)] from the dipole approximation to 
our general photon-matter interaction Hamiltonian, Eq. 
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Hed = ie^^^(nLjfe/2£oV^)'/'ekA • D„{6kAexp(ik • R) ~bl^exp{-ik ■ R)} \i) {j\ 

k A i,j 



(15) 



where \i) is the interacting charged particle state, or the eigenstate of Hm, R is the atom position, V = is the 
volume, T>ij = — e {i\ \j) are the matrix elements of atomic dipole moment, and A is the polarization. A general 
scattering transition rate r is 0|[Eq- (7.7.2)]: 



EE 



E 



{n - 1,1 J\ Hed \l) {l\ Hed \n,0,l) 



nto — oji 



27r 



(16) 



where |1) and |/) are the atomic ground state and final state. Initially, there are n photons with with frequency lj 
and wave vector k. At the end there are n — 1 incident photons and a single scattered photon with frequency lUsc and 
wave vector k^c- Then the cross-section follows from the relation cr(cj) = V/cnr, and the differential light-scattering 



cross-section is 



(in 



The differential cross-section of Rayleigh scattering is the special case when the atom returns 



to its ground state, which is ilJl[Eq. (8.8.1)]: 



dcr(w) 



4 4 

e uj 



E 



(e^c ■ Di;)(e • Dii) , (e • Di;)(esc ■ Dn) 



(17) 



where the parameters are obtained from matrix elements of Hed- 

The important equation is H16|l above: it shows that Rayleigh scattering originates from the interaction Hed- 
Clearly, the differential cross-section of Rayleigh scattering vanishes when Hed is zero. Our spatial BB method does 
just that: it effectively eliminates the interaction Hed- Of course, this is not unique to Rayleigh scattering, which is 
just one of the processes derived from considering various cases involving Hed- For example, photon absorption and 
emission is mainl y re lated to transitions involving two atomic or molecular levels. The corresponding matrix element 
for absorption is |ll][Eq. (4.10.1)]: 



"•kA - 1,2 



H 



ED 



n\^\, 1 ) = i%fcA exp[i(wo - Wfc)t + ik • R]n| 



1/2 
kA ' 



(18) 



where gi^\ = [eujk/2eQhV) ' - D12. The radiative lifetime is 

1/tr = 27r ^ ^ alxK^k 

k A 



^0) 



(19) 



and, of course, it follows from Eq. H18|l that this absorption is prevented when Hed is zero. 

Note how Hed, which is effectively eliminated by our method, involves the bosonic raising and lowering operators 
&kA, The reason that our method is so general is that it acts directly on these operators, and "time-reverses" 
Hed by flipping their sign. 



VIII. PROPOSAL FOR AN EXPERIMENT 



As mentioned above, a crucial requirement for the success of our proposed method is to insert the optical elements 
at intervals over which the coherence-loss is still quadratic, rather than exponential. An experiment to test for this 
regime is thus useful. This could be done by monitoring the coherence (in particular, loss) locally, by focusing onto 
the edge of the fiber and collecting light into a photon-counting device (since the absolute intensity would be very 
small). By moving the focus along the fiber, one should be able to track the decay as a function of distance from the 
fiber entry point, and observe the required quadratic-to-exponential transition, yielding an estimate of A. 

To actually test the method in the presence of phase-shifters, one could repeat the above experiment with a single 
fiber and write some phase-shift segments into it (as in the manufacturing of fiber Bragg gratings), at intervals 
bounded above by those determined from the first experiment. We note that a point of some potential concern is the 
impedance mismatch between air and the phase-shifter material, that will lead to reflection. Let {i = 1, 2) denote 
the indices of refraction: the reflected amplitude is (712 — ni)/(n2 -I- ni), which leads, at normal incidence, to 4% 
loss per air-glass interface. However, a standard anti-reflection coating can solve the problem: a quarter-wave layer 
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of material at ^rixn^ between the two materials (two equal reflections out of phase cancel out). In fibers the index 
changes will be smaller and reflection is typically neglected. Moreover, by writing a smooth phase profile as in the 
experiment proposed above, the reflection problem essentially disappears. 

Once A has been estimated, one can proceed to directly test our method, as follows. Take two fiber segments and 
write a tt phase-shifter (PS) into each. Attach them co-linearly (i) in the order PS-fiber-PS-fiber, (ii) in the order 
fiber-PS-PS-fiber, and perform a photon counting measurement. Our method should reduce attenuation in (i) by 
comparison to (ii). 



IX. CONCLUSIONS 



We have proposed a method to reduce quantum noise in optical fibers, via the insertion of phase-shifters at appropri- 
ately spaced intervals. We have shown that, in principle, this method can eliminate all quantum noise processes that 
do not involve photon number operators in the system-bath Hamiltonian; when such terms do arise, the phase shifters 
need to be supplemented with beam-splitters, and our conclusions remain. Thus, with simple linear-optical devices, 
quantum noise in optical fibers can be drastically reduced. This conclusion has potentially important implications for 
quantum communication (and its variants, quantum cryptography and distributed quantum computing) via optical 
fibers. The practical feasibility of our method hinges on the required distance between phase-shifters. We have given 
a rough upper-bound estimate of several meters based on known attenuation rates. We have also presented a more 
detailed calculation that predicts a range of distances, depending on the bath spectral density appropriate for a fiber. 
Ultimately we believe that the best way to test our proposal is to perform the relatively straightforward experiment 
that it implies. 
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APPENDIX A: DETAILED MODEL FOR ESTIMATING A 

Recall that our main approximation was the assumption of average fiber homogeneity, Eq. ((SJ. In this appendix 
we relax this assumption in order to estimate an upper bound on the distance A between phase shifters. We do this 
by considering corrections to order and the non-ideal case 

Ho{{k-l)A) = HoikA)+ePk, 

H\{{k-l)A) = H\{kA)+eQk, (Al) 

where e <C 1 and we take Pfc, Qk to be independent, identically distributed (IID) Gaussian, local, and time-dependent 
operator- valued corrections. This phenomenological model of fiber inhomogeneity may be the result of material non- 
uniformity along the fiber (such as local defects), slow time-dependent fluctuations in fiber properties, or even the 
quadratic interaction Q. By virtue of the central limit theorem it will be accurate in the case of a large number 
of defects. We assume that the effective BB time-interval r is chosen to be on the order of the small parameter e 
(though we make no attempt to estimate e). In this case, using the BCH formula e^e^ = e^+^+[^^B\/2+... second 
order (i.e., keeping only terms of order £, r, er, r^), we find instead of the ideal Eq. 



g-jH((fe-l)A)r-|-[g-iH(fcA)r-|-[ ^ g-iH((fe-l)A)r g-fflH(feA)nr 

^ g-j[Ho(feA) + H}(feA)+£(Pfc+Qfc)]Tg-i[Ho(feA)-H5(fcA)]r 

« e^^{-iT[2Ho{kA)+e{Pk + Qk)]-AH\{kA),Ho{kA)]}, (A2) 



where in the second line the effect of the phase shifters was to flip the sign (and thus cancel) the H\{kA) term. To 
the same order of accuracy the overall evolution operator becomes 

N/2 N/2 

U'{T,0) « e-^^«(°)^exp{-r2^[i/j(2A:A),i/o(2A:A)]}exp{-zer^(P2fe + Q2fc)}- (A3) 

fe=l k=l 

Let us evaluate the first exponential. Using Eqs. lO,(0: 

- i[H[{2kA),HQ{2kA)] = -i Y^[{Bl{2kA)bj + Bj{2kA)b]), nwj.{2kA){hj, + 1/2) + HM{2kA)] 

3,3' 

= ~i ^ {fujj{2kA)Bj{2kA) + [Bj{2kA), HMi2kA)]}b] 

3=0,1 

~{hwj{2kA)B]{2kA) - [B]{2kA),HM{2kA)]}bj 
= H', (A4) 

where H' is an effective Hamiltonian (it is Hermitian), which plays the role of a Lamb shift j2(]j|. We thus have for 
the first exponential in Eq. IjASp : 

N/2 

exp{-r2 Y,[Hi{2kA),Ho{2kA)]} = eyiv{-iT^H'), (A5) 
fe=i 

whose effect is an energy renormalization (i.e., a phase shift), and does not contribute to decoherence. 

Next, consider the second exponential in Eq. IjASp . The operator G{t) defined through X]^i(-f2fc + Q2k) ^ 
J^[P{t) + Q{t)]dt = G{t)dt is Gaussian distributed by our assumption that P2k,Q2k are Gaussian, IID random 
variables. We would like to estimate the average deviation in U'{T,0) that results from its presence. Since G{t) is 
Gaussian distributed the average can be computed as follows : 

N/2 



exp{-Z£T^(F2fc + Q2k)}) - (exp[— ier / G{t)dt)]) 
fc=i -^0 



exp[-ier / / {G{t)G{t'))dtdt'] 
Jo Jo 

exp[-eTr(r)]. (A6) 
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Expressed in terms of Fourier components of G{t) we have for the decoherence factor: 

2/ 



1 1"^ 

T{T) = - dLo{Gl)Q{co,T) (A7) 



where 

/2sin(.^/2)^^ (^3) 



Q{uj,T) = J^ dtdt' cos{uj{t-t')) ^ i^- 



UJ 



But in the Gaussian case we have (as in the spin-boson model pl|) 



(Gi)--/Hcoth^, (A9) 

where /(w) is the spectral density (of matter in the fiber) and (3 is the inverse temperature. Hence our result is that 
the correction is 



exp[-£Tr(r)] = exp 



-ET I dcji (cj) coth — 



/sin(wr/2)^ ^ 



UJ 



(AlO) 



The attenuation is thus strongly dependent upon the form of but also depends sensitively on temperature. In 

particular, the thermal time-scale hP is important in separating thermal effects from effects due purely to vacuum 
fluctuations |22| . In order to formally separate the two it is convenient to write 

coth ^ = 1 + n(w, /3) (All) 

where 

n(w,/3) = exp(-/3w/2)/sinh(/3tj/2) (A12) 

is the average number of field excitations at inverse temperature /3. 
In the limit of very low temperatures (/? ^ 1) we have 

n{uj, P) 2 exp(-/3cj) (A13) 

and we can analytically evaluate the integral in Eq. IjAlOp . e.g., for the class of Ohmic-type spectral densities, i.e., for 
the case 

/(w) = acj"e-"/"% (A14) 

where a is the coupling strength and tOc is the high-frequency cutoff (note that a is not dimensionless) . The result in 
the zero-temperature case is 



lim Jp°° dujl{uj) coth 



Puj fsm{ujT/2y ^ 



f log(l-f (cc>,r)2), n = l 

S.cj:?-ir(n - 1) (^1 - (1 + {uJcTf)"^ cos[(n - 1) arctan(cJcT)]) , n^l 



2 c 



To obtain the non-zero temperature correction in the approximation (|A13|I take these results, multiply by 2, replace 
ujc by i^p^^ everywhere, and add to the zero temperature case. We tabulate a few cases of interest in the zero 
temperature limit, letting x = lo^T: 

{(1 + x'^)-'^"!^, n = 1 (Ohmic) 

exp[-iae™cTf^], n = 2 (super-Ohmic) . (A16) 

exp[-2Q;e™^-p^^7j/J, n = 3 (Debye) 
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Let 1 — 5{T) be the desired coherence value after time T (or distance X); then we need to solve for the phase shifter 
spacing A from 

lim exp[-£Tr(r)] > 1 - S{T). (A17) 

We find (assuming a > 0): 

A2 < -4t>2ln[l-,5(T)]/ln[(l + a;2)], n 

The present model is, unfortunately, too phenomenological to make a reliable estimate of A. Nevertheless, it is of 
some interest to see its prediction. E.g., we could wish to improve upon the current figure of merit of 0.25 db/km to 
the threshold value of S(T) = 10"''. Recall that T = X/v, t = (A/w) and we assumed r ~ e. The coupling strength 
a is typically of order unity 0, H^; we shall set a = 1. We take v — c/1.6, the speed of light in a typical fiber, 
and 6{T) = 10~^. The results in the three cases, with x = 1.6/3 x 10~^Wc, are displayed in Fig.^ as a function of 
the high-frequency cutoff ujc- As a rough reference, the Debye temperature of amorphous silica is To = 342K [23 |. 
yielding a Debye frequency estimate of uic = kBTo/h = 2 x 10^'^Hz. The corresponding value of A is 6 x lO^m (n = 1), 
0.6m in — 2), lO^^m {n — 2>). This strong sensitivity to the decoherence model underscores the need for the proposed 
experiment in order to settle the question of the actual required distance between phase shifters. 



(A18) 



